ITP-SB-91-51 
Nov. 1991 



Completeness of the SO (4) extended 

Bethe Ansatz 
for the One-dimensional Hubbard Model 



Fabian H.L. Efiler 
Vladimir E. Korepin 
and 

Kareljan Schoutens 

Institute for Theoretical Physics 
State University of New York at Stony Brook 
Stony Brook, NY 11794-3840, U.S.A. 



Abstract 



We show how to construct a complete set of eigenstates of the hamiltonian of the 
one-dimensional Hubbard model on a lattice of even length L. This is done by using 
the nested Bethe Ansatz and the SO (4) symmetry of the model. We discuss in 
detail how the counting of independent eigenstates is carried out. 
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1 Introduction 



An important ingredient in the search for a theory of high-T c superconductors is the 
analysis of strongly correlated electron systems 0. A prototype model for these is 
the Hubbard model. Especially interesting is the //-pairing mechanism proposed by 
C.N.Yang in g g, |. 

The one-dimensional Hubbard model has been known to be exactly solvable 
since the work of E.Lieb and F.Y.Wu of 1968 6. In their paper, a large set of 
eigenfunctions of the hamiltonian were found by using the nested Bethe Ansatz @ . 
However, the issue of whether this set of eigenfunctions is actually complete has not 
been considered until recently. In a recent paper 0, we used the SO (4) symmetry 
of the Hubbard model (which has been explored in || [3|, f|, || |[) to show that the 
Bethe Ansatz is not complete (see below for some comments on this result). 

In this paper we will show that the 5*0(4) structure can be used to extend the set 
of Bethe Ansatz eigenstates to a complete set of eigenstates of the one-dimensional 
Hubbard model. This result was first announced in [10|. Here we provide a detailed 
account of the derivations involved. 



The Hubbard model describes electrons, which can hop along a one-dimensional 
lattice and which interact with coupling constant U if two of them occupy the same 
site. The hamiltonian of the Hubbard model on a periodic one- dimensional lattice 
of even, finite length L is given by (the small modifications in the potential term as 
compared to || are such that the 50(4) symmetry becomes explicit, see || |S], ^J) 

# = -E E (4,<rc i+ i,« + i^) + ^EKi - him,-! - ~) • (i-i) 

1=1,7=1,-1 i = l Z A 

Here the Cj jCr (cr takes the values ±1) are canonical Fermi operators on the lat- 
tice ii = 1, ...,L labels the lattice sites), with anti-commutation relations given by 
{c| i(7 , Cj tT } = Si t j5^ T . They act in a Fock space with the pseudo vacuum | 0) defined 
by Cj ;Cr 1 0) = 0. The operator n ita = c\ a Ci,(, is the number operator for electrons with 
spin a on site i. U is the coupling constant and can be either positive (repulsive 
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case) or negative (attractive case). For later convenience we define 

»=!■ ™ 

The analysis by E.Lieb and F.Y.Wu in || resulted in a large number of eigenstates 
of the hamiltonian, which are characterized by momenta fcj and rapidities A Q , where 
i = 1,2, . . . , N e and a = 1, 2, . . . , M for an eigenstate with a total number of N 
spin-up and M spin-down electrons. Our convention throughout the paper will be 

iV= number of spin-up electrons 

M— number of spin-down electrons 

N e = N + M= total number of electrons . 

We will now discuss the Bethe wave functions wave functions in the form as given 
by Woynarovich (Tl|, which is equivalent to the form found by Lieb and Wu. The 



nested Bethe Ansatz provides us with the following set of eigenstates with M spins 
down and N spins up : 

M N e 

\^ m ,n) = ^-i,„.,-m i( x i,---, x nJ H c i jt i n ckiio) 

l<x k <L j=l i=M+l 

= ^ CT1 , CT2 ,..., CTJVe (^n--- ) ^J IT 4,, CTj |o>, (1-3) 

l<x k <L j=l 

where we have put o\ = ... = om = — 1, &m+i — ■■■ — a N e — 1- 

The Bethe Ansatz wave functions explicitly depend on the relative ordering of 
the x.. We represent this dependence by a permutation Q of Af e elements, which 
is such that 1 < x„ < x r , < ... < x n < L. In the sector Q the general Bethe 
wavefunction for M spins down and spins up reads 

V'oi.oa,...,^ (^n-'-^jVe) =22 s 9n(Q) sgn(P) p i Q i <p(y 1: . . . , y M \P). (1.4) 

The P-summation extends over all permutations of N e elements and sgn(U) is the 
sign of the permutation II (II = Q,P). The amplitudes (f(y 1 , ■ ■ ■ ,y M \P) are of the 
form 

M 

cp( yi ,...,y M \P)= Y A,J{F P {K v yi) (1.5) 
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with 



/••/-(A,,//) = ( n ' 



where we defined 



and 



^ sin(k Pi ) - Aj -%\( 1 \ 



e2' ) (2) = ^(^)-A,±| (1.7) 



(1. 



By 7r = (7Ti, 7r 2 , .., 7r t , 7T t+1 , ..7Tm) we denote a permutation of M elements (spin-down 
electrons) and (t,t+i)7T = (^i, 7r 2 , .., 7r i+ i, tt*, ..7Ta/). A solution of ( |1.8| ) is given by 

^ n ( VV )- < L9 > 

\<Kk<M \ iv Tfc / 

The amplitudes ip(y 1 , ■ ■ ■ , V M \P) depend on o"i, a^ e and on Q through the numbers 
yi, ...,y m, which are defined to be the positions of the down spins among the spins 
in the series <TQ 1 ,aQ 2 , ■■■&Q Ne in increasing order, i.e., 

1 < yi < 2/2 < y 3 < <VM<N e . (1.10) 

For example, for one spin down and one spin up ( and a% = — 1, 02 = 1 ) we have 
the two cases y = 1 (if the spin down is to the left, which holds in the Q = (id) 
sector) and y = 2 (if the spin down is to the right, which holds for Q = (21)). 

As we already indicated in (|1.3| ), we will choose the notation such that the M 
down spins are at the positions x 1} ...,x u , i.e., a x = ... = a M = —1 and u M+1 = ... = 

° Ne = !■ 

We see that all solutions are characterised by N e momenta {kj\j=i,...,N e } of 
charged spinless excitations (holons), and M rapidities {Afc|jfe=i,...Af} of spin waves 
(spinons). 
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Imposing periodic boundary conditions on the Bethe Ansatz wavefunctions leads 
to the following equations for the parameters ki and A a 



ikjL TT \_JU « 4i_ 7 = 19 N 



*Z sm(Jb) -K-t ¥r A/3 - K - £ n , , 

sinffc ) — A + ^ = ~ Afl — A 4- — ' = U ' " " " M " 

i=l sm l K jJ /v « + 4i /9=1 A /3 ly <* + 2i 

Energy and momentum, i.e., the eigenvalues of the hamiltonian (|1.1|) and the 
logarithm of the translation operator, of the system in a state corresponding to a 
solution of ( |1.11| ) are 



E Ne = - 2 E C0S (^) + 17 
1=1 z 
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P = (1-12) 

8=1 

The second term in the expression for the energy in due to the shift of rij^ by | in 

Because of the antisymmetry of the product over c^'s under interchange of any 
two of them, the wave functions 4 , a 1 ...a Ne (x 1 ...x N J can be (and have been) chosen 
to be completely antisymmetric under the simultaneous exchange Xk < — > Xj and 
°k < — ► 0j, i-e., 

' l PcT 1 ..a : j..a k ..a Ne (x i ---X j ...X k ...X Ne ) = — 1p ai ,.a k ..aj..a Ne \ X i ■ ■ - x k ■ ■ - x j ■ ■ - x Ne ) • (1.13) 

We now define 'regular' Bethe Ansatz states (for finite L), to be denoted by 
I ^m,n)i by the properties that N — M > (non-negative third component of the 
spin), N e < L (less than or equal to half filling), and that all A a and all kj are finite. 
Bethe Ansatz states with N — M < and/or N e > L can be obtained from the 
regular Bethe Ansatz states by using simple symmetry operations, i.e., reflection of 
the third component of the spin and particle/hole correspondence, which commute 
with the hamiltonian @. The model is invariant under spin- rotations, with the 
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corresponding SU (2) generators given by 



C = £<&Vi' C ] = (0\ c* = ~X>,-i-^i)- (i.i4) 

i=l Z 8=1 

(Note that ( z equals minus the third component of the total spin.) For even L the 
model has a second SU(2) invariance, generated by 

L L 1 L L 

^ = 0-1)^,1^,-1, ^ = XX-i)H-i4i> ^ = ;EK-i+»i,i)- «■ ( L15 ) 
i=i i=i z i=i z 

The raising operator ry of this second SU(2) creates a pair of two opposite-spin 
electrons on the same site, with momentum tt. Combining the two 5f/(2)'s, which 
commute with the hamiltonian and with one another, leads to an 50(4) invariance 
of the one-dimensional Hubbard model for even lattice leng thsll. For a discussion 
of the theoretical and possible experimental consequences of the existence of this 
symmetry, which also exists in the Hubbard model in two or three dimensions, we 
refer to the papers @, [|, J], [13| . 



In a previous paper 0, we established the following remarkable property of the 
regular Bethe Ansatz eigenstates of the hamiltonian: they are all lowest weight states 
of the 50(4) algebra (pip , flEIgP , i.e., 



V I ^Pm,n) =0, CI ^m,n) = 0- (1.16) 

This shows that acting with the raising operators rf and on | iPm,n) leads to 
new eigenstates of the hamiltonian, which are not in the regular Bethe Ansatz. In 
this way, every regular Bethe Ansatz state | iPm,n) is the lowest weight state in a 
multiplet of states, which form a representation of 50(4). Since 

1 1 
Vz | V>M,iv) = ~{N e - L) | iPm,n) , Cz | iPm,n) = -(M - N) I ^ m ,n) , (1.17) 

a state \^m,n) nas s pi n ?7 = \{L ~ N e ) with respect to the ?7-pairing SU{2) algebra 
and spin £ = ^(N — M) with respect to the ( SU(2) algebra. The dimension of the 
corresponding 50(4) multiplet is therefore given by 

dim M)N = (2»7 + l)(2C+l) = (L- N e + l)(N - M + 1) . (1.18) 
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The states in this multiplet are of the form 

By symmetry, all the states that are highest or lowest weight states with respect to 
one of the SU(2) algebras are again given by the Bethe Ansatz (although in general 
they are outside the regular Bethe Ansatz). All other states are not given by the 
Bethe Ansatz, which shows that for this model the Bethe Ansatz is not complete. 
The simplest example of a state that is outside the Bethe Ansatz is | 0), which 
describes a single repair of momentum tt. The fact that this state is outside the 
Bethe Ansatz was explicitly confirmed in 0. 

It is the main purpose of this paper to show that, if one counts the number 
of eigenstates that are related to the regular Bethe Ansatz states by the SO (4) 
symmetry, one finds 4 L , which is precisely the correct dimension of the Hilbert 
space of the model 0. Thus we will conclude that the Bethe Ansatz together with 
the 5*0(4) structure leads to a complete set of eigenstates of the one- dimensional 
Hubbard model. 

The paper is organized as follows. In section 2 we discuss in some detail the 
nature of the solutions of the Bethe equations for the Hubbard model. In Appendix 
A we show that the so-called A and k — A strings give wave functions that describe 
bound states. In section 3 we will then count eigenstates and prove completeness. 
An explicit construction of the |(L + 2)(L — 1) solutions of the Bethe equations in 
the sector with one spin-up and one spin-down electron (as opposed to the somewhat 
indirect construction used in the general proof) is presented in Appendix B. 

2 Solutions of the Bethe equations for the Hub- 
bard model 

Let us focus on the Bethe Equations (|1.11[) , which express the fact that the Bethe 
Ansatz wavefunctions ( |1.4j ) satisfy periodic boundary conditions. 

*There are 4 possible electron configurations per lattice site (spin up, spin down, both spin 
up and spin down, and empty site), thus the corresponding direct product Hilbert space is 4 L 
dimensional. 
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Counting regular Bethe Ansatz states means counting inequivalent solutions of 
the equations fll.ll ) while taking into account the 'regularity conditions' N — M > 



and N e < L. Following M.Takahashi o, we will first distinguish different types of 
solutions {ki,A a } of ( |1 . 1 1|) . The idea is that for a solution {ki,A a }, the set of all 
the fcj's and A a 's can be split into (three) different kinds of subsets ('strings'), which 
are 

1. a single real momentum fcj 

2. m Aq,'s combine into a string-type configuration ('A-strings'); this includes the 
case m — 1, which is just a single real A a < 

3. 2m fcj's and m A^'s combine into a different string- type configuration ('/c-A- 
strings') |. 

For large lattices (1 << L), almost all the string configurations are close to 'ideal- 
ized' string-solutions where both the /c's and the A's are assigned imaginary parts 
according to a 'equal-spacing' prescription @. For a A-string of length m the ra- 
pidities involved are 

K' j =K-\(rn+l-2j)u Aureal j = l,2,...,m. (2.1) 



The /c's and A's involved in a /c-A-string are 





= 71 — 


sm -1 (A 


kl 


= sin" 




kl 


= 71 — 


k 2 


kl 


= sin" 


-l (Ar _ 


kl 


= 71 — 


7.1 

ot 



1 , 

-mu) 

2 ' 



-(to-4)w) 



^these correspond to bound states of spin waves (magnons) 
Hhe case m = 1 describes a 'Cooper pair' of electrons 



and 



k 2 a m ~ 2 = sm- 1 (A' a m + ^(m-2)u) 

k 2 a m = ir-wT 1 (A!™+~mu). (2.2) 

A^ = A7_I( m + l_2j>, A'J" real j = l,2 m. (2.3) 

]3D are valid up to exponential corrections of order 0(exp— 5L), where S 
is real and positive (and depends on the specific string under consideration). In 
Appendix A we discuss the wave- functions corresponding to some of the the string 
configurations 2. and 3. and show that they correspond to bound states. 

Let us now consider a solution that splits into M m A-strings of length m, M' n 
fc-A-strings of length n (containing 2n k^s and n A a 's) and M e additional single fcj's. 
Clearly, we have 

oo oo 

N e = M e + 2 £milC, M= ]Tm(M m + ilC). (2.4) 

m=l m=l 

How many solutions of this type exist? 

The idea is that each of the strings in a solution can be characterised by the 
position of its center (a real number), which we denote as in ( |2.1| ) by A™, a = 
1,2,... M m , for the length-m A-strings, by A'™, a = 1, 2, . . . M' m for the length-m 
fc-A-strings (as in ( |2.3| ) ) and which is simply equal to kj for the unpaired momenta 
kj, j = 1, 2, ... , M e . Because of the periodic boundary conditions, these parameters 
will have to be chosen from a discrete set. 

Following [fll] , we now write the following equations for the parameters kj, A^ 1 
and A^™. They follow from ( p. . 1 1| ) and the form of the 'idealized' string-solutions 
which we discussed above (we write M' = J2m=i m M' m ) 

oo M n 

k j L = 2nI j -J2Y, d 

n=l a=l 

N e -2M' / ah • i, \ oo M m / A n A m ^ 
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sin 1 ( A'™ + m^j + sin 1 f A£ - in j 



N e -2M' 



'A'™ — sin fcj 
n[7 



oo M^, / A ' n 

771=1 /3=1 V 



A /3 



C7 



, (2.5) 



where 



6{x) 



2tarT 1 (4x) 



4 x 



n+m— 2 



4 x 



n+m, 

for n m 
for n = m 



(2.6) 

The Jj, J™ and J^ n are integer or half-odd-integer according to the following pre- 
scriptions: Ij is integer (half-odd-integer) if X) m (M m +M^) is even (odd); the J™ are 
integer (half-odd-integer) if (N e — M n ) is odd (even); the J'™ are integer (half-odd- 



integer) if [L — (N e — M' n )) is odd (even). According to fl4j| , we have the following 
inequalities 

1 oo 

1^1 < -(N e -2M'-J2tnmM m -l), 



m=l 



J':\ < (L-N e + 2M'- Y,tnmM' m -l) 



0</ J < 



m=l 



(2.7) 



where t nm = 2 Min(n, m) — 5 nm . 

We will now make the standard assumption that, in order to enumerate the 
different solutions of the system ( |2.5| ), it is sufficient to enumerate all possible sets 
of non-repeating (half-odd) integers Ij, J™ and J'™, satisfying ( |2.7|) . 

[This assumption mimics the similar assumption which is usually made for the 

It is known, however^, that the actual 



L Heisenberg XXX model 



16] 



spin- 2 

distribution of the different types of solutions can be different from the one implied 



§This fact was actually already noticed in the original paper by Bethc, fVt 
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by this counting. We carefully studied this phenomenon in a recent paper fl8 
where we give a detailed discussion of the two-magnon sector of the XXX model. 
We explicitly show that the deviations from the above assumption can be viewed as a 
'redistribution phenomenon', which does not affect the total number of Bethe Ansatz 
states. For the Hubbard model we have a similar situation, which gets complicated 
further by the fact that there is a free coupling constant U in the model. In Appendix 
B we analyze in detail the N = 1, M = 1 sector of the Hubbard model. Although 
we do find ^/-dependent redistributions among different types of solutions, we find 
agreement with the predictions based on the 'ideal' string assumption for the total 
number of states in this sector, which is \{L — 1)(L + 2) for the L-site model.] 

From ( p. 71) we read off that the numbers of allowed values for the (half-odd- 
) integers corresponding to each of the fundamental strings are 

1. L for a free ki 

2. N e - 2M' - Em=i t nm M m for a A-string of length n 

3. L - N e + 2M' - E^=i t nm M' m for a k- A-string of length n. 

The total number of ways to choose the (half-odd-) integers in a solution with mul- 
tiplicities M e , M m and M' m is therefore given by (remember that the integers are 
assumed to be non-repeating) 

, ^ , AT e - 2M' - E^-i t nm M^ , 
n{M e ,{M m },{M' m })=\ ITT , x 




M, 



n 




Ne+ 2M' -E: =1 t„ m i: 

ML 



(2- 



n 



The total number of solutions of ( |1 . 1 1| ) with given numbers N and M is now 
obtained by summing n(M e , {M m }, {M' m }) over all the M e , M m and M' m , under the 
constraints (|2.4|). 

Every solution to (p. .11 ) gives us a regular Bethe Ansatz state, which comes with 



an entire multiplet of eigenstates of the hamiltonian, the dimension dmiM.iv of which 
is given in ( |1.18|) . The full number of eigenstates that are obtained from the Bethe 
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Ansatz and the SO (4) symmetry is therefore given by 
# (eigenstates) = 

i- oo oo oo 

EE E E E n(M e ,{M m },{M^}) 

M>0N>0 1 M e =0M m =0M^=0 
N-M>0 N+M=M e +2 V°° , mM' 

N+M<L N-2M=M e -2Y°° , mM m 

— L — ^ m — 1 



dim M . 



N 



(2.9) 



The counting of the eigenstates that are obtained from the SO (4) extended 
nested Bethe Ansatz has thus been reduced to a purely algebraic problem, which 
we will solve in the next section. 



3 Counting eigenstates 

In this section we will prove that for general even L the sum in (|2.9| ) equals 4 L . This 
will prove completeness. Before we come to that, we show the examples of the 2-site 
and 4-site models. The 2-site model (L = 2) was discussed in |7j , where we presented 
the explicit form of a complete set of 4 2 = 16 eigenstates of the hamiltonian. In 
Table 1 we show how the counting presented in section 3 works out in this case. 



M e 






M 


N 


n 


dim Mj7V 


# (states) 

















1 


3 


3 


1 











1 


2 


4 


8 


2 











2 


1 


3 


3 


2 


1 





1 


1 


1 


1 


1 








1 


1 


1 


1 


1 


1 










16 



Table 1: L = 2. n denotes the number of regular Bethe Ansatz states of a given 
type. There are a total number of 16 = 4 2 eigenstates of the hamiltonian. 



The total number of 16 states splits into 2 singlets, 2 triplets and 2 quadruplets 
of SO (4). The ground state is the singlet with Mi = 1 for the case U > and the 
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Mi 


M 2 


Mi 

1 


ML 


M 


N 


71 


dim M jv 


# (states) 
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5 


5 


1 

















1 


4 


8 


32 


2 

















2 


6 


9 


54 


3 

















3 


4 


8 


32 


4 

















4 


1 


5 


5 


2 


1 
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1 


i 

1 


6 


3 


18 


U 


U 


U 


1 


U 


1 


1 


Q 
O 


Q 
O 


n 

y 


3 


1 











1 


2 


o 
o 


A 


oz 


1 








1 





1 


2 


8 


4 


32 


4 


2 











2 


2 


1 


1 


1 


4 





1 








2 


2 


1 


1 


1 


2 


1 





1 





2 


2 


6 


1 


6 











2 





2 


2 


1 


1 


1 














1 


2 


2 


1 


1 


1 


4 


1 











1 


3 


3 


3 


9 


2 








1 





1 


3 


6 


3 


18 










256 



Table 2: L — 4. There are 60 regular Bethe Ansatz states, which, when weighted 
with the correct 5*0(4) multiplicities, give a total of 256 = 4 4 eigenstates of the 
hamiltonian. 
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singlet with M[ = 1 for the case U < 0. In both cases it is a bound state of one 
spin up and one spin down electron with energy Eq = — Ju 2 /4 + 16. The counting 
for the 4-site model (L = 4) is presented in Table 2, where we show how the total 
number of 4 4 = 256 is obtained. (Notice that the total number of regular Bethe 
Ansatz states is only 60.) 

We now turn to the proof that for general (even) L the sum in ( [2.9| ) equals 4 L . 
We will split this proof into two steps as follows. In the first step we will prove the 
following two identities 



£ n 

Mi,Af 2 ,...=0 n= \ 
Y2=i mM ™= M 





(3.1) 



and 





For later convenience we define 



(N -2M+1) = 2 



N 



m=l 



(3.2) 



(3.3) 



and 



oo / P n 
n({M m }) = II 

n=l \ M n 



(3.4) 



where t nm = 2 Min(n, m) — 5 nm as before. In the second step we will then use 
identities (|3.1|) and ( |3.2|) to perform the summation in ( |2.9| ). 

The auxiliary identities (|3.1|) and ( p.2|) have a natural interpretation in the con- 
text of the spin- 1/2 Heisenberg XXX model jlB]] . The equation ( |3.1|) gives the total 
number of regular Bethe Ansatz states (defined by M < [iV/2]) with M overturned 
spins in the XXX model on a lattice of length N. The second formula shows that 
the total number of states obtained by combining the regular Bethe Ansatz with 
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the SU(2) structure equals 2 N , which is the dimension of the Hilbert space of the 
XXX model. These relations thus establish the completeness of the SU{2) extended 
Bethe Ansatz for the XXX model. 

The fact that identities that have their origin in the XXX model play a role 
here should not come as a surprise. Indeed, our method of solution of the Hubbard 
model is the nested Bethe Ansatz. The solutions to the Bethe Ansatz are specified 
by two sets {kj} and {A a } of spectral parameters. The k/s are momenta of charge 
density waves, whereas the A a 's, which describe the 'nesting' of the Bethe Ansatz, 
are rapidities of spin density waves of the type encountered in the Heisenberg XXX 
model. This should make clear that the second stage of the nested Bethe Ansatz 
for the Hubbard model is really a spin-problem, which is very similar to the Bethe 
Ansatz analysis of the Heisenberg XXX model. Our two-step procedure for per- 
forming the summation is natural from the point of view of the nesting: in the first 
step we sum over the spin degrees of freedom, and in the second step we then sum 
over the charge degrees of freedom as well. 



STEP 1. 



Let us now explain how the equation (|3.1|) can be derived. In the first step, one 
simply solves for Mi = M — J2m=2 m M m and substitutes this back into the left hand 
side of ( [3.1|) . Using this value for M 1; the quantities P n reduce to 



Pt = N - £ t lm M m = N - M +J^(m- 2)M m , 

m=l m=3 

oo oo 

P n = N-J2 tnmMm = N — 2M + M n + 2 ]T ( m - n ) M ™ ■ (3-5) 



m=l m=n+l 



Let us now consider the summation over Mi. Although our summand in the left 
hand side of ( [3.1|) has the form of an infinite product, only two of the factors contain 
the variable Mi. Singling these out, one finds that the summation over Mi is as 
follows 



n 2 = £ 

M 2 =0 



/ N-2M + M 2 + 2 E^ =3 (m - 2)M„ 
v Mi 
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N-M + E™ =3 (m-2)M„ 

x I I . (3.6) 

In order to perform this summation we will make use of the identity 

E( B a a )^ = (i-*r x - fl . ( 3 - 7 ) 

which can easily be proved by induction. As a simple consequence, we have 



wr'-a+aou- ), (3.8) 



0=0 



where the notation | x a in the left hand side means that we single out the coefficient 
of the power x A . The right hand side of (|3.6j ) is of the same form as (|3.8|) and we 
find 

n 2 = (l_^ ) -l-[Ar-2M + 2(M3+2M 4+ ...)] (l+x) iV-M + (M3+2M 4+ ...) 

00 

= (1 + x) n - m {i - x 2 )-^ 1 ' 1 n (4 o) ) M - u 

= ^ / ^T^ 1 + ^) iY " A/ ( 1 - x 2 )-™^ 1 ft (^°>)^ , (3-9) 



where 



(1 _ X )2(n-2)( 1+x )n-2 " ^ LU ) 



In the last line of (|3.9|) we extracted the coefficient at x M by performing a contour 
integral around the origin x — 0. After performing the M2 summation ( |3.1| ) now 
reads 



Mi,M 2 ,...=0 
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where 

A(x) = (l + xf^il-x 2 )- 1 '™ x 



oo 



/ N - 2M + M n + 2 E~ =n+1 (™ - n)M^ 

En „ n(^ (0) ) Mi -(3-i2) 



oo oo 

X 

M 3 ,M 4 ,...=0 n=3 



V 



The summation over M 3 is given by 



=3 




N - 2M + M 3 + 2 E^ =4 (n - 3)M n \ 
Q 3 = E I ^ J 

= (L^)- 1 "^- 1 ^. (3.13) 
At this point, the full expression for -A(x) has been reduced to 

-1-N+2M 



A(x) = (1 + x)^- M (l - x 2 )- 1 -™ (l - Z f ] ) + x 
( N — 2M + M n + 2E» =n+1 (m - n)M„ 



oo oo 

x e n 

M 4 ,M 5 ,...=0n=4 



V 



(3.14) 



where i^ 1 ) are defined through 



4 1 ) = " . (3.15) 

(i-zf) 2(n " 3) 

From the above it is now clear, that the sum with respect to M 4 and all M n with 
n > 4 has the same structure as the sum with respect to M 3 . Thus the final result, 
after performing all summations, will look like 

A(x) = (1 + x) JV - M F(x)- 1 - iV+2M , (3.16) 

where 

oo 

F(x) = (l-x 2 ) I] (l-^ 3) ) (3-17) 

m=3 



17 



and we have the iteration formula 



57(m-l) 

4 m) = ~ Ti 9T • (3-18) 

An+2 ) 

Our task is now to find a closed expression for F(x) by exploiting this relation. 
We define 

U 2 = x~ 2 , U m = -±^, m>3, (3.19) 

Zm 

so that F(x) can be written as 

F W = fi (l - jrl ■ (3-20) 

We now claim that the functions U m (x) satisfy the following recursion relation, to 
be denoted by RR I p 

RRI P : (U p+ z — l) 2 = U p+ aU p+ 2 , p>0. (3.21) 

Together with the initial conditions 

£,„ = ,-», C f,= ' 1 -^ 1 + ^ (3.22) 

these relations completely fix the functions U m (x) and thereby the function F(x). 

In order to prove the recursion relation ( |3.21| ), we first give a second recursion 
relation, which involves some of the other Z's and which we shall denote by RR II p 

RRII P : ^gi = ^, p >0, n>p + 3. (3.23) 

Zk ' U P+3 

Let us now show that the validity of both recursion relations can be proved by 
induction. We start at the point where we have Z^\ which is defined by ( |3.10|) , 
and U2 and U3 as above in ( 3.22Q . One easily checks that RR II p= o is valid. Using 



( gig ) for the definition of U A = l/z[ 1] and RR II p =0 , one proves RR I p =q. This 
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establishes the validity of both RR I p=0 and RR II p= o, which is the starting point 
for the proof by induction. 

Let us now assume that we have proved the validity of both RR I p and RR II p 



for some given p. By using this induction assumption and the definition ( 3.18Q , one 
then proves the relation RR II p+ i ( |3.23| ). After that, by using the definition ( |3.18| ) 
and RR IIp+%, one then proves RR I p+ \. This completes the induction step. We 
may thus conclude that the relations RR I p and RR II P are valid for all p > 0. 
One easily checks that the expressions 

U^M^-jfS 1 ) 2 (3-24) 



a(x) — a(x) 



with 




satisfy the recursion relations (|3.21|) and the initial conditions Q3.22j ). The function 
F 2 (x) is now expressed as a convergent product 



f(x)= n— ^— = n— ^ — = }™w 2 ar = a(x) x(x + lh (3 - 26) 

m=2 m m=2 m ^ * 



where we defined Ui(x) = in accord with Q3.21| ), and where we used (p. 21 



m 



the second equality. This brings us to the following representation of the number 
of regular Bethe Ansatz states with M overturned spins (using (|3.1|), ( |3.11| ), ( |3.12| ), 
([U6p, ( gg§ , ( |3T25D and (|3~2p ) 



E 

M m =0 
M=V°° . mM„ 



n{{M m }) 



1 r dx 



,N-M 



'(1 + X) + yJ(l+x)(l-3xY 



2M-N-1 



2m J y M+1 



2dy (2(l + y)f- M 



1 + y+Jl-y 



2M-N-1 



(3.27) 
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where the contour is a small circle around the origin and we used the substitution 
y = tz~ ■ Calling y^ 1 = cosh0, the integral reduces to I + — J_, where 

'* = hr + = (/ M ' Tl ) (3.28) 

where we wrote = A — iip, with A — > oo. This finally establishes the result ( |3.1|) 
since 



E 



n({M m }) 



M m =0 

L^dra — 1 



AT-1 \ _ / iV— 1 
.V -M 1 I V -• M 1 



M J 1 M-l 



(3.29) 



We still have to prove (|3.2|), which can be done as follows 



[N/2] 

E 

M=0 



N 
M 

[N/2] 



N 
M-l 



(N-2M+1) 



W*\ / >r \ [iV/2]-l / v 

E m ("-2M+1)- E m (^-2M-1) 

M=0 \ / M=0 \ / 



[JV/2]-l 

2 E 

M=0 
[JV/2]-l / 

E 

Af=0 \ 
JV 

E 

M=0 



N 
M 



+ 



N 

N 
2 



N 
M 



N 

N 
2 



N 
M 



l + iV-2 



AT 



l+JV-2 



N 
~2 



N 

+ E 

JV-[iV/2-l] 



A^ 
M 



(3.30) 



This completes the proof of equation ( |3.2j ). 



STEP 2. 
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The total number of states that are obtained from the SO (4) extended Bethe 
Ansatz for the Hubbard model is given by ( |2.8| ) and (|2.9| ). The summations over the 
multiplicities M m and over the difference N—M in the summation ( j2.9| ) are precisely 
of the type fl3TT|) and fl3~2j), respectively, if we substitute M — > \[M e — N + M) and 
N — > M e . (Under these summations the total number of electrons N e is kept fixed.) 
The summation that remains after this 'spin summation' is 



# (eigenstates) 



E (L~N e + l) 

N e =0 



X 



N e oo 

E E : 

M e =0 ML=0 



L 

M P 



n 

71=1 



L- N e + Y^ = i^m-t nm )M' m 



M' 

n 



(3.31) 



where as before t nm = 2 Min(n, m)—5 nm . In the next step we perform the summation 
with respect to all M^'s, using a similar kind of 'summation device' as in STEP 1. 



As a consequence of (|3.7|) we have 



x 2 r l - B (\ 



2 X y 




2P x 2M[+p 



(3.32) 



and therefore 



1 
2~n% 



dx 



x 



7+1 



-1-B 



(l + 2xY 




2M[ 



2^~ 2M i. (3.33) 



The integration is along a small contour around zero. Defining E = L — N e , 7 = 
-2YZ,=2 mM 'm and B = E + 2£^ =2 (m - the r.h.s. of flQ3[) becomes 



the summation over M[ in ( 3.31 ), if we solve the constraint in the sum in (|3.31|) 
for M e = N e - 2£~ =1 mM' m . Using (|3~33| ) in flOID we then obtain the following 
expression for the number of eigenstates : 



# (eigenstates) = / 

2ixi J 



L+1 f* 2 , (1 + 2*) L E (E + 1)7^^ 

.T^ +i (l — X 2 ) ^ (1 — X z )^ 



F(x), (3.34) 
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where 



Fix) 



oo oo 

e n 

m>2 



( E + M' n + 2Y.™ =n+ M-n)M' m v 

n K 0) 

m=2 



(3.35) 



and 



2m 



(1 - X 2)2(m-1) 



The summations over M%, M^,... have precisely the form of the l.h.s. of 
can thus be performed easily. The result is 



(3.36) 
and 



F{x) = j[{i-zsr*>r l - B 



m=2 



(3.37) 



where 



zip- 1 ) 



(i - z^f^-p- 1 ) ' 



It can be shown along the lines given in STEP 1, that the quantities U, 
obey the recursion relation 



(3.38) 



i 



m _(m-2) 



(Up+2 - l) 2 = C/p+sC/p+i , V > o 



(3.39) 



with initial conditions 
Equation (|3.23|) is replaced by 



crip) TT 

_ u p+l 

criP) ~ TJ 



'1 — X 



2\2 



X 1 



p > 0, n > p + 2 . 



(3.40) 



(3.41) 
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Equation ( |3.34|) now can be written as 



1 r rlx — 
# (eigenstates) = — j ^(1 + 2x) L £ (E + l)x E [f{x)]- E - 1 , (3.42) 



dx 

-{l + 2xr 

E=0 

where 

oo 

f(x) = U^-Ur 1 ). (3.43) 
i=i 



The solution of the recursion relation ( p. 39 ) is again of the form ( p. 24 ), i.e., Uj 

( a(xy+ 1 -a(x)-i- 1 \ 2 i 

( L{x)-a(x)-i ) > Wllere nOW 



<*) = Tx + ^- 1 (3 ' 44) 



due to the new initial conditions (|3.40|) . Insertion of the resulting expression for U\ 



into ( p. 43 ) leads to the following result for the function f(x): 



2f(x) = 1 + VI - 4x 2 = 2x a(x) . (3.45) 
Equation ( |3.42|) can now be rewritten as 

L 

# (eigenstates) = + 1) , (3.46) 

_B=0 

where 

^ = 2s/ d B) (l + ^I'Wr-'- (3-47) 

The contour integration can be worked out as in section 2. Defining a = A — itp 
with A ^> 1 and substituting x = ea + e - a we obtain 

1(E) = I + (E) - I_(E) , (3.48) 

where 

J± (£) = J_ e ±« e -(i+*)« ( e f +e ~f) 2L . (3.49) 

27T JO V ' 
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Expanding 



(ef+e-f) 2L = E( 2 p L )e^-*> (3.50) 



=o 



and then using 

±- ^ dtp e*™ = 5 n>0 (3.51) 

Z7T JO 

in the resulting expression, we find that 

'^-(l-e)- '-W = (l-b- 2 )- < 352 » 

Plugging these results into (|3.48f) and then ( |3.46|) we are left with only a single 
summation 

# (eigenstates) = £ (J3 + 1) j ^ L — E ) _ ( L — E — 2 ) } ' (3 ' 53) 



This summation can be performed the same way as ( |3.30|) and we finally obtain the 
desired result 

# (eigenstates) = 4 L . (3.54) 



This concludes our two-step evaluation of the sum Q2.9Q 



Using the above, we can obtain a closed expression for the number of regular 
Bethe Ansatz states for given numbers M and iV of spin-down and spin-up electrons: 



oo oo oo 



E E E n{M e ,{M m },{M' m }) = 

Af e =0 M m =0 M' m =§ 
N e =M e +2j2 mM L 

(a^)((m) + (m-2))~((a^ + i) + (a^-i))(m-i)' 

(3.55) 

This formula is the close analogue of the result (|3.1| ) for the XXX Heisenberg model. 
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We repeat once more our conclusion, which is that the combination of the nested 
Bethe Ansatz with the SO (A) symmetry of the one-dimensional Hubbard model leads 
to a complete set of 4 L independent eigenstates. 

It is a pleasure to thank C.N. Yang for proposing the ideas worked out in this 
paper. We thank S. Dasmahapatra for stimulating discussions. This work was 
supported by NSF Grant PHY- 9107261. 
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A Bound states in the one-dimensional Hubbard 
model 



In this appendix we investigate the nature of the A and k — A strings in the one- 
dimensional Hubbard model. We show for explicit examples, that both kinds of 
strings lead to certain kinds of bound states (i.e., the wave function decays expo- 
nentially with repect to the differences of coordinates) . 

A.l A strings 

Let us consider the example of N electrons with spin up and two electrons with spin 
down forming a A string, i.e., 

A 2 = Ai* , ki, real . (A.l) 

The set {Ai, A 2 | fci, kN+2} must fulfill the periodic boundary conditions (|1.11|) . 
The wave function corresponding to this set of spectral parameters is given by 

ip_ 1 _ 111 (x 1 ,x 2 ,...,x N+2 )=Y,sgn(Q)sgn(P)e 1 -^ p ^ XQ ^(y l ,y 2 \P) , (A.2) 

Pes 2+N 

with amplitudes 

ip( yi ,y 2 \P) = A id F P (Ai,yi) F P (A 2 ,y 2 ) + A 21 F P (A 2 , yi ) F P (A 1 ,y 2 ) 



A e£>(Pi) / e M(P n ) \ ?i e?>(Pi) ) e?(P, 



in • 



!J2 



- 1 e«(P*) \ 1 / KJZ 1 e ®(Pj \ 1 



n ^ "TnTTTT n 



in. - 



(A.3) 

We want to show that this wave function decays exponentially with respect to 
the difference of the coordinates y\ and y 2 . The only nontrivial y dependent part of 



27 



the wave function are the amplitudes <p. Therefore it is sufficient to prove that they 
decay exponentially. 

Taking exponentially small corrections A into account, the A string is of the 
form 

u + A 



Ai = A + 



. . u + A . 
A 2 = A W ith A* = -A 



(A.4) 



As all momenta k* are real this leads to the following inequalities 



< 1 . 











>1 , 




e«(i) 





(A.5) 



Using the periodic boundary conditions 



A 2 - Ai - u _ Ad 
A 2 - Ai + u ~ A 2X 



(A.6) 



we can express A 2 ± in terms of A id . Using the second set of periodic boundary 
conditions 



e W (j) eL 2) (j) 



(2), 



(A.7) 



we can express products over e± in terms of exponential factors of magnitude 1. 
Straightforward computations yield 



v ( yi ,y 2 \P) = A ld e^ k -i L [ 1] 



V2-1 ( •-' (/, 



r-ui. e +\Pi) j e +\Pyi) e +\Py2) 



+ 



yi - 1 e®(P t ) 



1 



1=1 



i=V2 e + (Pi) J e + (^2) e +\Pyi, 



(A.8) 



By definition y 2 > y\ and the inequalities ( A.5 ) ensure that the factor in brackets 
in the first term (and thus the whole term) in (|A.8|) decays exponentially for y 2 ^> y%. 
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The second term can be dropped, because for spin waves N plays the role of the 
lattice length, and in order to investigate asymptotic properties of the wave function 
we should consider the infinite volume limit, i.e., N — > oo. In this limit the second 

term can be set to zero as ( n£j/2 Va)!^! ) vanishes. 
A. 2 k — A strings 

We consider the example N = M = 2, i.e., 2 electrons with spin up and 2 with spin 
down. The periodic boundary conditions read 

e »,L _ £Hj)^U) 4 

A fc — A string solution of these equations takes the following form in the L — ► oo 
limit: 

Ai = a-| = a; 

ki = 7r — arcsin(A — u) 
k 2 = arcsin(A) = n — k 3 

k 4 = 7r — arcsin(A + u) . (A. 10) 

In the finite volume there exist two distinct configurations that both lead to k — A 
strings in the limit L — ► oo, depending on whether ki and k 3 are real or complex 
for finite L. 

Case (i): k 2 and k 3 are complex 

In this case the fc's can be rearranged such that Im(ki) < Im{k-2) < < 
Im{k 3 ) < /m(/c 4 ) and the invariance of the periodic boundary conditions under 
complex conjugation gives the additional constraints 

k* = k 2 , k* A = h . {All) 
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Taking this into account one obtains the following solution of ( |A.9| ) for finite but 
large L: 

a, = a-^ = a; 



sin(ki) = A — u — — + €i = sin(kl) 
sin(k 2 ) = A + - + e 2 = sin(k^) , 



(A.12) 



where S is purely imaginary and £^2 are complex. The exponentially small corrections 
are of the orders 



5 = O (2u e -^ +k ^ L ) , e x = Q (-2u e~ ifclL ) , e 2 = O (-2m e~ ifclL ) . (A.13) 



The wave function is given by ( |A.2j ) with N = 2 and ( A.3 ). After inserting the values 
of the spectral parameters found in (|A.12 ) into ( |A.2|) and (|A.3| ) and re-normalising 
the resulting expression one finds : 

sgn(Q)(-l)yi + y* (e^U^i - j^U^A if (y u y 2 ) £ {(i j2 ), (3,4)} 



else 



(A.14) 



where R is the permutation (1, 3, 2, 4). 

Due to the ordering of the imaginary parts of the momenta kj both terms on the 
l.h.s. of ( |A.14|) decay exponentially with respect to the magnitudes of differences of 
coordinates \xp. — Xj\ and thus the wave function describes a bound state. 

Case (ii): k 2 and k% are real 



In this case we must drop the constraint k\ 
solution to ([A.9Q 

u + 5 



k 3 . We then find the following 



Ai = A 



A* 
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sin(ki) = A— u — - + e 1 = sin(k* A ) 
sin(k 2 ) = A + e 2 



sin(k 3 ) 



A + e 3 



(A.15) 



where e-y is complex and e 2j 3 are real, while 5 is again purely imaginary. The correc- 
tions are of the orders 



The computation of the wave function is analogous to case (i), the only difference 
being a new re- normalisation constant. The wave function is given by the same 
expression as in case (i). Again it describes a bound state although k 2 and k 3 are 
now real. 

B The M = N = 1 sector in the Hubbard model 

In this appendix we further work out the stucture of the Bethe Ansatz wavefunctions 
in the sector M = N = 1. In that sector the wavefunctions depend on parameters k±, 
k 2 and A. Our general analysis in sections 3 and 4 gives the following possibilities: (i) 
we can have M e = 2, Mi = 1, which gives real values for k±, k 2 and A, or (ii) we can 
have M[ = 1, which gives a m = 1 /c-A-string with complex conjugate ki, k 2 and real 
A. According to the counting of section 3, we expect to have \L{L — 1) real solutions 
(i) and L — 1 string solutions (ii), giving a total number of \(L — 1)(L + 2) Bethe 
Ansatz states in this sector. However, we already mentioned that, in the context of 
the XXX Heisenberg model, there is a 'redistribution phenomenon' between different 
types of solutions, which does not affect the total number of states within a sector, 




5 






(A.16) 
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but which does affect the distribution of those states over various types of string- 
solutions and real solutions @. In this appendix we will establish a very similar 
result for the Hubbard model: we will find that the numbers of the real solutions 
(i) and string solutions (ii) are not always given by the values quoted above, but 
that the expected total number \{L — 1)(L + 2) of solutions in this sector can be 
rigorously established. 

Let us consider the Bethe equations for periodic boundary conditions in the 
sector M = N = 1 

^k lL = sin fci — A — | ^ L = sin k 2 - A - f 
sin h - A + | ' sin k 2 - A + f 

Since we assume A to be finite, we can solve for it and find A = ~(sin&i + sin/^). 
The equations then reduce to 

eikiL = sinfci-sinfc 2 -tt ^ e i {kl+ k 2 )L = l _ / B 2 x 

sin ki — sin k 2 + u 

We can solve the second equation by putting k\+ k 2 = ^m, with m — 0, 1, . . . , 2L— 1 
and write k\ = |Jm + x and k 2 = jm — x. The remaining equation reads 

p i(n m+Lx) = -^cos(^)sinx-i 

-^cos(^)sinx + r 1 J 

One easily checks that, if Xq solves this equation for m = m , then x = Xq + n 
solves the equation for m = mo + L, and that the resulting wavefunctions are the 
same. We can thus restrict our attention to m = 0, 1, ... L — 1 and —it < x < tt. 

Let us now try to find real solutions x for the equation ( p.3| ) for given m. Taking 
a logarithm we have 

/ U .irm. \ 1 , n 
arctan I — — cos(— —) sinil = --Lx — im , (,h>.4J 

where n is an arbitrary integer for m odd and half an odd integer for m even. It is 
rather straightforward to solve this equations by a graphical method: one plots both 
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the l.h.s and the r.h.s. (for various n) of these equations on the interval — tt < x < tt 
and reads of the intersection points, which are then solutions of the equation. For 
large enough \U\ this procedure is easily carried out and one finds the following. 

For m even there are solutions for n = — . . . , which are L solutions in 
total. Since solutions x and — x are equivalent (and x = is not among the solutions) 
we should divide this number by 2. Using that there are L/2 possible even values 
for m we thus find L 2 /4 solutions. For m odd one finds non-equivalent solutions 
for n = 1, 2, ... , for each m, which gives a total number of (L 2 — 2L)/A. (The 
solutions x = 0, x = ±7r, which exist for generic U, give vanishing wavefunctions 
in general.) Adding up the contributions from odd and even m, we find \L{L — 1), 
which is indeed the number predicted by the counting in section 3. 

However, let us now assume that m is odd and that U is close to a critical value 
U m , which we define by 

U m ,nm L 

— cos(— ) = -. (B.5) 

At the value U = U m the curve for the l.h.s. of ( p. 4]) has slope L/2 at x = ±7r, 
and at the value U = —U m the curve for the l.h.s. of ( |B.4| ) has slope L/2 at x — 0. 
Since the r.h.s. is given by straight lines of slope L/2, and since both curves already 
had intersections at x = (for n = 0) and x = ±7r (for n = ±L/2) (which did 
however not give rise to non-trivial wavefunctions), it will be clear that the number 
of intersections changes when U reaches the critical values ±[/ m . In fact, one finds 
one extra real solution x (together with the equivalent solution —x) for a given odd 
m as soon as U < \U m \. For example, if U is such that U\ > U > U3 >....> there 
will be one extra real solution to the equations (|B.2|). 



The complex values for x which solve the equation ( |B.3|) are of the form x = iy or 
x = n+iy with y real. In a way similar to what we showed above, one can analyse the 
equations for the real quantity y by a graphical method. If \U\ is sufficiently large, 
one finds precisely one complex solution for m = 1,2, ... ,L with the exception of 
m = L/2. In that case, there are thus L — l complex solutions, which is in agreement 
with the counting of section 3. 

However, from the graphical analysis one finds that the complex solution for a 
given odd m disappears as soon as \U\ is chosen to be smaller than U m . Note that 
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this happens precisely in the regime where we have found one extra real solution! 

We thus find a redistribution phenomenon, where solutions change their nature 
as a function of U, in analogy to what we found for the XXX Heisenberg model 
in | 18| . In the Hubbard model the phenomenon is easily understood: if \U\ is made 



small enough, the interactions become so weak that some of the bound states (with 
complex x) decay into real solutions (with real x). 

When \U\ is chosen to be equal to one of the critical values U m , there do exist 
nontrivial wavefunctions with x = or x = ±7r, i.e., with coinciding k\ and k 2 . These 
wavefunctions can be seen to be nonvanishing by a renormalisation a la l'Hopital. 

In all cases, the total number of eigenf unctions of the hamiltonian in the sector 
M = N = 1 is found to be \{L — 1)(L + 2), which is the value predicted by the 
counting in section 2, and used for the proof of completeness in section 3. 
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